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ABSTRACT 

We report the detection of a new 3.06 ms binary pulsar in the globular cluster NGC6544 using a 
Fourier-domain "acceleration" search. With an implied companion mass of ~ 0.01 M Q and an orbital 
period of only Pf, ~ 1.7hours, it displays very similar orbital properties to many pulsars which are 
eclipsed by their companion winds. The orbital period is the second shortest of known binary pulsars 
after 47 Tuc R. The measured flux density of 1.3 ± 0.4 mJy at 1332 MHz indicates that the pulsar is 
almost certainly the known steep-spectrum point source near the core of NGC6544. 

Subject headings: galaxy: globular clusters: individual: NGC6544 — pulsars: individual: 
PSR J1807-2459 — radio continuum: stars 



1. INTRODUCTION 

Globular clusters are rich sources of radio pulsars. Since 
the discovery of the first pulsar in a cluster by Lyne et al. 
(1987), a series of deep searches have accumulated almost 
50 pulsars in globular clusters, with the majority of these 
being binary millisecond pulsars (MSPs). See Lyne (1995) 
and Kulkarni & Anderson (1996) for reviews. 

Fruchter & Goss (1990; 2000, hereafter FG00) used the 
Very Large Array (VLA) to detect steep-spectrum radio 

. sources, which are likely pulsars, in numerous globular 
clusters. They discovered a relatively bright point-like 
source within 5" of the optical center of cluster NGC6544 
(1.2±0.1mJy at 20cm and 6.5±0.8mJy at 90cm). Biggs 
& Lyne (1996) searched NGC6544 at 50 cm as part of a 

. survey of globular clusters and found no pulsars down to 
a limiting flux density of ~ 4 mJy. 

It has been argued that Doppler smearing from orbital 
motion could explain the lack of detection of pulsations 

' from the pulsar-like sources found in VLA surveys (FG00). 
The sensitivities of observations whose durations exceed 

■ even a few percent of the orbital period may be drasti- 
cally reduced for traditional analyses (Johnston & Kulka- 
rni 1991). Various schemes have been developed to cor- 
rect, at least partially, for this loss of sensitivity when 
the Doppler effect of the orbital motion can be approxi- 
mated as a constant frequency derivative (e.g. Middleditch 
& Kristian 1984; Anderson et al. 1990; Wood et al. 1991). 
These "acceleration" searches have met with significant 
success, including the discovery of nine new binary MSPs 
in the globular cluster 47 Tucanae (Camilo et al. 2000). 



We have discovered a binary MSP in NGC6544 using 
Fourier-domain techniques to correct for pulsar orbital mo- 
tion. This pulsar is almost certainly the radio source re- 
ported by FG00 near the center of the cluster. Results 
from our single observation indicate that this object shares 
many of the properties displayed by the burgeoning class 
of eclipsing binary MSPs. These systems are characterised 
by short orbital periods (P& ~ 1 — 10 hours) and very low- 
mass companions (m2sini ^ 0.05 M©) (Rasio et al. 2000; 
Nice et al. 2000). D'Amico et al. (2000) independently 
discovered this pulsar and report the results of follow-up 
observations. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed NGC6544, a dense core-collapsed globu- 
lar cluster, on 1998 February 7 with the multibeam data 
acquisition system on the Parkes radio telescope in Aus- 
tralia. Signals from each of two orthogonal linear polar- 
izations were measured and summed from 256 contiguous 
0.25 MHz wide channels centered at 1332 MHz. Each chan- 
nel was one-bit sampled at 8 kHz and written to magnetic 
tape. The full observation comprised 28.9 minutes of data. 
Staveley-Smith et al. (1996) and Lyne et al. (2000) discuss 
the observing system in detail. 

Since no pulsars were known in this cluster before our 
observation, we searched a wide range of possible disper- 
sion measures (DMs) for candidate signals. We dedis- 
persed the data into time series of 13865600 points us- 
ing 600 trial DMs from to 600pccm~ 3 , in increments of 
lpccm~ 3 . This stepsize allowed a worst-case dispersion 
smearing over our bandwidth of ~ 0.1ms. After barycen- 
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tering the data using the JPL DE200 ephemeris (Standish 
1982), we performed a Fast Fourier Transform (FFT) on 
each time series and saved the resultant spectra to disk. 

Traditional time-domain acceleration searches are per- 
formed by stretching or compressing a time series to 
compensate for a constant frequency derivative and then 
Fourier transforming the resulting series. We have devel- 
oped a Fourier-domain acceleration search (S. Ransom & 
S. Eikenberry 2000, in preparation) to compute the co- 
herent Fourier response over portions of the frequency- 
frequency derivative (/— /) plane using only local Fourier 
amplitudes from an FFT of the whole dataset. Similar 
techniques were described in Middleditch et al. (1993). 

The method correlates a predicted Fourier response with 
a subset of the complex Fourier amplitudes in the initial 
full-length FFT. We effectively apply a matched filter to 
the data which "sweeps up" all of the signal power that 
the orbital motion spreads over nearby frequency bins. 

Fourier-domain acceleration searches offer several signif- 
icant advantages over their time-domain counterparts: 

1. Fourier methods do not require stretching or 
compressing the time series. Time-domain 
stretching or compressing is usually performed 
using linear interpolation, which changes the 
statistics of the data. Portions of the data are 
effectively two-bin averaged which decreases 
sensitivity to high frequency signals. 

2. Time-domain techniques require a full-length 
FFT for each trial acceleration. For very long 
integrations (which are common in cluster 
searches) that do not fit into the computer's core 
memory, FFTs take orders-of-magnitude more 
time to compute than usual (e.g. Bailey 1990). The 
Fourier technique requires only a single full-length 
FFT for each observation (or trial DM). 

3. Correlations using only localized Fourier am- 
plitudes are fast and memory efficient. The 
correlations can always be performed in core 
memory and are calculated efficiently using short 
FFTs and pre-computed response templates. The 
memory locality also allows efficient parallelization 
of Fourier-domain searches. 

4. The Fourier-domain method allows the cal- 
culatation of only independent /-/ trials. 
With time-domain techniques, as described by 
Camilo et al. (2000), any choice of acceleration 
stepsize results in either under- or over-sampling 
accelerations for the vast majority of frequencies 
searched. This amounts to loss of sensitivity or 
wasted CPU cycles respectively. A small amount 
of over-sampling can be incorporated to eliminate 
"scalloping" in both the / and / directions (see 
e.g. van der Klis 1989). 

We performed Fourier-domain acceleration searches on 
each of the long FFTs looking for signals that drifted by 
up to 100 Fourier frequency bins during the observation. 
If we define z to be the number of Fourier bins that the 
pulsar frequency (/ Q ) drifts over the course of an integra- 
tion, the corresponding average acceleration of the pulsar 



is a = zcj '{T 2 /„), where T is the total integration time. 
For this observation, a frequency drift of 100 Fourier bins 
corresponds to an acceleration of 500 ms~ 2 for a 10 Hz 
signal or 5 ms -2 for a 1000 Hz signal. Since most known 
binary pulsars show maximal accelerations of only a few 
ms~ 2 , we were sensitive to all but the most exotic binaries 
or pulsars with spin periods (P sp in) much less than ~ 2 ms. 
Pulse profiles of strong candidates were folded at a series of 
trial DMs, frequencies, and frequency derivatives around 
the acceleration search values in order to maximize the 
signal-to-noise. 
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Figure 1 — An 18 a (single trial) detection of the 1.7 hour binary 
PSR J1807-2459 in the globular cluster NGC6544 using a Fourier- 
domain "acceleration" search. Contour intervals correspond to 30, 
60, 90, 120, and 150 times the average local power level. The intrin- 
sic pulsar period and / = (which corresponds to an un-accelerated 
FFT of the data) are marked by the solid black lines. The dashed 
line is the predicted "path" of the pulsar given the orbital solution 
in Table 1. During the 28.9 minute observation, the pulsar moved 
from ~ 11 o'clock to ~ 3 o'clock on the ellipse. The peak's slight 
offset from the ellipse and the presence of "shoulders" indicate that 
the constant / assumption of the acceleration search could not fully 
correct for the orbital motion during this observation. 



3. RESULTS 

We detected a strong pulsar candidate (power/average 
local power ~ 165 corresponding to a single trial signifi- 
cance of 18 a) at a DM of 134pccm~ 3 . The average fre- 
quency was 326.85895Hz and the signal had drifted by 8.6 
Fourier bins (a = 2.6 ms~ 2 ) during the observation (see 
Figure 1). Search techniques that do not account for ac- 
celeration would have detected the signal at about ~ 34 
times the average local power level, corresponding to an 
8 a detection. 

Optimization of the candidate by pulse-folding produced 
a narrow pulsar-like profile at a localized maximum in DM 
of 134pccm~ 3 indicating that the signal was almost cer- 
tainly not due to terrestrial interference. The Taylor & 
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Cordes (1993) free-electron model gives an estimated dis- 
tance of 3.3 kpc (~ 25% error) for this DM in the direction 
of NGC6544, in fair agreement with the published cluster 
distance of <~ 2.6 kpc obtained by fitting the mean V mag- 
nitudes of horizontal branch stars (Harris 1996). 

We estimated a pulsar flux density at 1332 MHz of 
1.3 ± 0.4 mJy by comparing the integrated pulsed flux 
and observed noise to the predicted system noise level 
given a system sensitivity and total system temperature 
°f T tot ~ T sys + T s k y . The multibeam system has a cold- 
sky T sys ~ 21 K and a sensitivity of 1.36 JyK -1 . Haslam 
et al. (1982) measured a sky temperature at 408 MHz 
in the direction of NGC6544 of T sky ~ 260 K. Assum- 
ing a spectral index of —2.6, which is typical for this re- 
gion (Lawson ct al. 1987), we obtain T s k y ~ 12 K and 
Ttot ~ 33 K at 1332 MHz. FG00 reported a 20 cm flux 
density of 1.2 ± O.lmJy and a spectral index of —1.1 be- 
tween 90 cm and 20 cm which corresponds to a flux density 
of <~ 1.4 ± O.lmJy at 1332MHz. Our measured flux den- 
sity agrees with that of FG00 within the errors, suggesting 
that the pulsar may be the previously unidentified point 
source near the core of NGC6544. Since the DM-based 
distance estimate for the pulsar matches published values 
for NGC6544, cluster membership is very likely. 

The relatively large signal acceleration indicates that the 
pulsar is in a short period binary with a low-mass com- 
panion. These systems, almost without exception, have 
eccentricities very near to zero due to tidal circularization 
(see e.g. Phinney & Kulkarni 1994 and Rasio & Heggie 
1995). 
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Figure 2 — Orbital fit for PSR J1807-2459 assuming an ec- 
centricity of zero. The error-bars plotted are ±2 a uncertain- 
ties in the pulse arrival times. An eclipse, if present, should 
occur near True Anomaly = 90°, which corresponds to su- 
perior conjunction. The apparent lack of an eclipse could 
be due to the relatively high observing frequency (1332 MHz). 
The inset plot shows the observed pulse profile after fold- 
ing with the best-fit orbital ephemeris and pulsar frequency. 
The error-bar indicates ±2<r uncertainties in the profile values. 
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Figure 3 — Error ellipses (from inner to outer of 1 <r, 2 tj, and 3 a) 
for the orbital period and semi-major axis of PSR J1807— 2459. The 
large covariancc is due to the fact that the observation covered only 
~ 30% of the orbital period. 

In order to study the assumption that the pulsar is a 
member of a low-mass binary with an eccentricity of zero, 
we split the observation into 16 equal-length parts and 
folded the data in each at the nominal pulsar frequency. 
We determined phase offsets for each of the pulse arrival 
times and performed a Levenberg-Marquardt least-squares 
fit of a circular orbit (i.e. sinusoid) to the phase offsets. 
The results are tabulated in Table 1 and the residuals and 
error ellipses of the fit are shown in Figures 2 and 3 re- 
spectively. 



Table 1 




Parameters for PSR J1807-2459 


Parameter 


Value 


Right Ascension a (J2000) 


18 h 07 m 20?36(2) 


Declination a (J2000) 


-24° 59' 52'/6(4) 


Dispersion Measure (pc cm~ 3 ) 


134(2) 


Flux Density at 1332 MHz (mJy) 


1.3(4) 


Pulse FWHM at 1332 MHz (%) 


11(3) 


Pulsar Period (s) 


0.003059447(2) 


Epoch (MJD) 


50851.9 


Orbital Period (days) 


0.070(3) 


Projected Semi-Major Axis (lt-s) 


0.0116(8) 


Eccentricity 13 


0.0 


Epoch of Periastron (MJD) 


50851.887(1) 


Mass Function (M ) 


3.4(5) x 10~ 7 


Companion Mass Limit c (M©) 


> 0.009 



a A. Fruchtcr 2000, private communication. 

Fruchter & Goss discovered an error in the 20 cm 
position reported in FG00. The corrected position 
shown above will appear in an erratum. 



b Parameter assumed to be exactly zero. See §3. 
c Assuming a pulsar mass of mi = 1.4 M Q . 
Note — Numbers in parentheses represent 2 a 
uncertainties in the last digit quoted. 
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We found no other candidates in our search of NGC6544 
to a limiting flux density of ~ 1.1 mJy for 1 ms period sig- 
nals and ~ 0.5 mJy for periods > 10 ms. 

4. DISCUSSION 

There are currently eight binary pulsar systems 
that display similar orbital characteristics to those of 
PSR J1807-2459: PSR B1957+20 (Fruchter et al. 1988) 
and PSR J2051-0827 (Stappers et al. 1996) in the Galac- 
tic disk, PSR J1910+0004in NGC6760 (Deich et al. 1993), 
and pulsars 47 Tuc I, J, O, P, and R (Camilo et al. 2000). 
Each of these systems displays a pulse period of a few 
milliseconds, an orbital period of 1-10 hours, and a very 
low-mass companion (m 2 sini ;$ 0.05 M©). Some of these 
systems display radio eclipses, particularly at longer wave- 
lengths (i.e. 50-90 cm). 

Our current 1332-MHz data show no evidence for an 
eclipse (see Figure 4) at or near superior conjunction, but 
we cannot rule out the possiblity of eclipses at longer wave- 
lengths. Pulsars J2051— 0827 and 47 Tuc J show eclipses at 
radio wavelengths between 50-90 cm, but are detected at 
all orbital phases at 20 cm, scintillation permitting (Stap- 
pers et al. 1996; Camilo et al. 2000). PSR J1910+0004 
in NGC6760 was detected at all orbital phases at 20 cm, 
but Deich et al. (1993) could not rule out short duration 
eclipses or eclipses at longer wavelengths. 

Most formation theories for short-period binary MSPs 
in globular clusters involve dynamical interactions with 
primordial binaries which have been shown to exist in sig- 
nificant numbers (Hut et al. 1992). Rasio et al. (2000) 
describe a recent variant of the standard "spin-up" model 
(see Phinney & Kulkarni 1994, for a review) using an ex- 
change interaction between a neutron star (NS) and a hard 
primordial binary containing at least one relatively mas- 
sive main sequence (MS) star (~ 1-3 M©). If the exchange 
interaction results in a binary containing the NS and the 
massive MS star, the system enters a common-envelope 
(CE) phase once the MS star evolves and fills its Roche 
lobe. At the end of the CE phase, the NS emerges with a 
very low-mass companion in a very short period circular 
orbit. Systems with orbital periods Pb 8 hours undergo 
evolution due to gravitational radiation and enter a second 
phase of mass transfer that spins the NS up to millisecond 
periods. 



The scenario described by Rasio et al. (2000) is possi- 
ble only if the exchange interactions occurred at a time 
when massive MS stars were still present in the cluster. 
The NSs and primordial binaries must also have undergone 
mass segregation and ended up near the cluster core before 
the MS stars evolved. The timescale for mass segregation 
is the half-mass relaxation time, which for NGC6544 is 
~ 200 Myr (Harris 1996). This is significantly less than 
the MS lifetime of a 1-3 M© star and implies that this sce- 
nario may have produced PSR J1807— 2459. 

Alternatives to "spin-up" models exist in which the NS 
is created either by accretion induced collapse of a white 
dwarf (WD) (e.g. Bailyn & Grindlay 1990), or the coales- 
cence of a pair of massive WDs (Chen & Leonard 1993). 
Pairs of primordial binaries interact to produce hard bi- 
naries which contain a massive WD or WDs that these 
models require. The NSs that result from these models are 
born spinning rapidly (P sp in ~ 2 — 10 ms), have relatively 
weak (10 9 — 10 10 G) magnetic fields characteristic of MSPs, 
and may have low-mass companions in very short period 
orbits. Interestingly, these systems could be the progeni- 
tors of the low-mass X-ray binaries (LMXBs) rather than 
their descendants (Chen et al. 1993). 

Future observations of PSR J1807-2459 will lead to im- 
proved measurements of rotational and orbital parameters, 
an estimate of or upper limit to the orbital eccentricity, 
and perhaps information about the circumstcllar medium 
if eclipses are detected. Rasio et al. (2000) predict that 
many similar systems must exist in the globular cluster 
system — some with orbital periods as short as ~ 15min. 
If recent successes such as the 20 cm searches of 47 Tu- 
canae (Camilo et al. 2000) are any indication, this pre- 
diction will soon be put to the test. The ever-increasing 
speed of computers and improved algorithms for binary 
pulsar detection (e.g. Ransom 2000) will allow the analy- 
sis of extremely long observations (i.e. days or weeks) that 
should reveal even the weakest of pulsars in the tightest of 
binaries. 

We would like to thank Vicky Kaspi, Josh Grindlay, Gio- 
vanni Fazio, and Nichi D'Amico for useful discussions. FC 
is supported by NASA grant NAG 5-3229. SE is supported 
by a NSF CAREER grant. 
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